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Abstract. A new instability mechanism is described in 
accretion flows where the gas is accelerated from a station- 
ary shock to a sonic surface. The instability is based on a 
cycle of acoustic and entropic waves in this subsonic region 
of the flow. When advected adiabatically inward, entropy 
perturbations trigger acoustic waves propagating outward. 
If a shock is present at the outer boundary, acoustic waves 
reaching the shock produce new entropy perturbations, 
thus creating an entropic-acoustic cycle between the shock 
and the sonic surface. The interplay of acoustic and en- 
tropy perturbations is estimated analytically using a sim- 
plified model based on the compact nozzle approximation. 
According to this model, the entropic-acoustic cycle is un- 
stable if the sound speed at the sonic surface significantly 
exceeds the sound speed immediately after the shock. The 
growth rate scales like the inverse of the advection time 
from the outer shock to the sonic point. The frequency 
of the most unstable perturbations is comparable to the 
refraction cutoff, defined as the frequency below which 
acoustic waves propagating inward are significantly re- 
fracted outward. This generic mechanism should occur in 
Bondi-Hoyle-Lyttleton accretion, and also in shocked ac- 
cretion discs. 

Key virords: Accretion, accretion disks - Hydrodynamics 
- Instabilities - Shock waves - Binaries: close - X-rays: 
stars 



1. Introduction 

Hydrodynamic instabilities in the process of accretion 
onto a compact star may help to understand the time 
variability of the emission from X-ray binaries. For exam- 
ple, the supersonic motion of a compact star in a uniform 
gas, first described by Hoyle & Lyttleton (1939), Bondi & 
Hoyle (1944), was later found to be unstable in numerical 
simulations (Matsuda et al. 1987, Fryxell & Taam 1988, 
Matsuda et al. 1992, Ruffert & Arnett 1994). In such flows. 
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the gas captured by the accretor is first heated and decel- 
erated to subsonic velocities through a shock, and then 
accelerated to supersonic velocities towards the accretor. 
Observationnal consequences of the instability were dis- 
cussed by Taam et al. (1988), Livio (1992), De Kool & 
Anzer (1993). Nevertheless, this instability is still poorly 
understood even in its simplest formulation. A recent at- 
tempt by Foglizzo & Ruffert (1999) has led us to look in 
more detail at the interplay of acoustic and entropy waves 
in the subsonic flow between the bow shock and the accre- 
tor. We identify in the present paper an entropic-acoustic 
cycle: entropy perturbations advected with the flow trig- 
ger acoustic waves propagating back to the shock, where 
in turn they trigger new entropy perturbations. The linear 
stability of this cycle depends on the efficiencies of these 
two processes, measuring: 

(i) the amplitude of the pressure perturbations Sp~ 
propagating against the stream, triggered by the advection 
of an entropy wave of amplitude 5 Si , 

(ii) the amplitude of the entropy perturbation 682 trig- 
gered by the interaction of an acoustic wave Sp~ with a 
stationary shock. 

The entropic-acoustic cycle is a priori unstable if 682 > 
SSi. As will be seen throughout this paper, a complete 
description of the entropic-acoustic instability requires to 
take into account acoustic waves Sp'^ propagating with the 
stream, triggered by the interaction of Sp'' with the shock 
and refracted out by the flow gradients. 
Since such entropic-acoustic cycles may be present in as- 
trophysical flows other than the Bondi-Hoyle-Lyttleton 
accretion, the present paper is dedicated to a description 
of this generic mechanism in its most general formula- 
tion. Similar entropic-acoustic cycles have been studied 
extensively in the context of jet nozzles, where the grav- 
itational potential plays no role. In particular, the first 
process (i) is well known among the combustion commu- 
nity since the works of Candel (1972) and Marble (1973). 
A clear overview of the subject can be found in Marble & 
Candel (1977). An example of entropic-acoustic instabil- 
ity in nozzles is the ramjet 'rumble' instability studied by 
Abouseif, Keklak & Toong (1984), where the combustion 
zone plays the same role as the shock in the second pro- 
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Fig. 1. Illustration of the release of acoustic energy by the 
advection of a localized perturbation of entropy in a nozzle 



cess (n). 

The paper is organized as follows: in Sect. we analyze 
the excitation of acoustic waves by the advection of en- 
tropy perturbations. The entropy perturbations produced 
by a sound wave reaching a shock are estimated in Sect. 0. 
These two processes are combined into a global mode in 
Sect. 0, where we characterize its spectrum of eigenfre- 
quencies. Possible astrophysical applications are discussed 
in Sect. 0. 



2. Acoustic energy released by the advection of 
entropy perturbations in an irrotational flow 

2.1. Energy released by the advection of a localized 
perturbation of entropy 

We wish to estimate the energy 6E12 released by the ad- 
vection of a perturbation SS of entropy into a gravitational 
potential, from a region I where the unperturbed velocity, 
sound speed and gravitational potential are vi,ci,^i to 
a region II closer to the accretor where their respective 
values are V2,C2,^2- These quantities are related through 
the conservation of the Bernoulli constant B in the un- 
perturbed flow. Denoting by 7 the adiabatic index of the 
gas, and by h its enthalpy, 



7-1 



(1) 

(2) 



A direct calculation of SE12 is possible when the unper- 
turbed flow is uniform before and after the region of con- 
vergence, so that the entropy perturbation has enough 
time to reach a pressure equilibrium with the surround- 
ing gas (see Fig. |^). Let m -I- 6m be the mass of the gas 
with a perturbed entropy S + SS. From the conservation 
of entropy. 
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For the sake of simplicity, the ratio of the molecular weight 
^ to the gas constant TZ is set to /i/7?. = 1 throughout 
this paper, with no loss of generality. The total energy 
5Ei2 released during the advection through the nozzle is 
calculated in Appendix A: 



5Ei2 = (/i2 - hi)5m. 



(4) 



This calculation, which involves no linearization, indicates 
that the the release of energy during the advection of a 
localized entropy perturbation is a fundamental process 
associated to the variation of enthalpy in the flow. 5Ei2 
corresponds to the sum of the energies of two localized 
pressure perturbations propagating with and against the 
stream. According to Eq. (Q), the bigger the increase of 
entropy through the flow, the bigger the release of acoustic 
energy by the advection of entropy perturbations. In an 
accretion flow governed by a gravitational potential <& = 
GM/r, the lengthscale of the inhomogeneity of the flow 
scales like the distance r to the accretor. Consequently 
there is no region II where the flow is uniform as in Fig. nl 
Nevertheless, the time needed by an entropy perturbation 
to reach a pressure equilibrium with the surrounding gas 
is of the order of the sound crossing time 51/ c, where 51 
is the length of the entropy perturbation along the flow. 
The distance it travels before the pressure equilibrium is 
reached is of the order of M51, where M. = v/c is the 
mach number of the flow. The above calculation of 5Ei2 
is thus relevant for an accretion flow in the region where 
the length of the perturbation is short enough compared 
to the scale of the flow inhomogeneities: 



A^^«l. 

r 



(5) 



Note that Eq. (^ does not specify the relative ampli- 
tudes of the pressure perturbations propagating with and 
against the stream. 

2.2. Acoustic flux released by the advection of an entropy 
wave 

If the entropy perturbation were an extended entropy 
wave, we expect acoustic waves to be triggered by the 
enthalpy gradient, by the same process as described in 



Sect. 2.1, By definition, the acoustic energy of an acous- 



tic wave is the part of the energy which depends (in a 
quadratic way) only on the linear approximation of this 
wave (e.g. LighthiU 1978, Chap. 1.3). The flux F± of this 
acoustic energy is defined by: 



F^ 



Moc- 
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(6) 



where Mq is the mass flux and A1 = v/c the Mach number. 
The index (-I-) is used to denote the wave propagating in 
the direction of the flow, (— ) otherwise. Any perturbation 
can be decomposed in the linear approximation onto the 



T. Foglizzo & M. Tagger: Entropic-Acoustic instability 



acoustic and entropic modes as a sum of a pressure pertur- 
bation Sp = Sp'^ + Sp~ and an entropy perturbation SS 
(e.g. Landau & Lifshitz 1987, Chap. 82). The hnear re- 
sponse of the flow to small perturbations is described by 
the two complex coefficients TZi, Qi defined at a position 
ri in the flow as follows 

Pi Pi 

The coefficient TZi is directly related to the refraction of 
the acoustic flux of an incoming pressure perturbation Spf 
by the flow inhomogcneities. Using Eqs. (o) and (R) with 
SS = 0, the fraction of refracted acoustic flux is: 



1-Mi 



Ff {l + Mi 



TZi 



If there is no incident acoustic wave {5pi = 0) , the acous- 
tic flux Fj~ propagating against the stream is directly re- 
lated to the amplitude of the incident entropy perturba- 
tion, and depend on the cocfflcicnt Qi as follows: 



p- - aV 2 (1--Mi)^ ^2 (5S^ 



Both TZi and Qi depend on the frequency uj of the per- 
turbation. An accurate calculation of 7?,i, Qi in a realistic 
accretion flow is beyond the scope of the present paper, 
since it would require solving the full set of linearized Eu- 
ler equations. This will be done in the particular case of 
a radial flow in a forthcoming paper (Foglizzo 2000). The 
frequency dependence of TZi can be anticipated by remark- 
ing that the refraction is negligible for acoustic wave with 
a wavelength much shorter than the scale of the inho- 
mogcneities of the flow. In an adiabatic accretion flow, 
the gravitational potential $ = —GM/r is responsible for 
both the acceleration and the heating of the flow. The 
lengthscale of the flow inhomogcneities decreases towards 
the accretor, so that we expect a strong refraction at low 
frequency and a negligible refraction at high frequency. If 
the accretion flow is accelerated up to supersonic veloci- 
ties through a sonic surface, acoustic waves with a high 
enough frequency to penetrate in the supersonic region 
cannot be refracted out: a frequency cut-off uJcxa therefore 
corresponds to the highest frequency above which refrac- 
tion becomes negligible. 

2.3. Estimates ofTZi, Qi based on the compact nozzle 
approximation 

We build in Appendix B a toy model to obtain rough 
estimates of the entropic-acoustic coupling in an accretion 
flow where the mach number increases from A^ i < 1 to 1 . 
The refraction point ra{oj) of acoustic waves of frequency 
to separates the flow into two regions: 

- a region of propagation r > Va where both the 
entropic-acoustic coupling and the acoustic refraction are 
negligible. 



- a region r < ra where entropy perturbations are cou- 
pled to acoustic waves, and where acoustic waves are re- 
fracted. 

The propagation of acoustic waves in the flrst region is 
easily described using the conservation of their acoustic 
flux. The region of coupling is described using the "com- 
pact nozzle approximation" , used for jet nozzles by Marble 
(7) & Candel (1977). This approximation is valid for pertur- 
bations with a long enough wavelength in order to treat 
the nozzle as a discontinuity in the flow. The formulae 
obtained by Marble & Candel (1977) are extended in Ap- 
pendix B to the case of a flow in a potential $. In this aca- 
demic case, the nozzle separates the upstream region char- 
acterized by a uniform velocity Va, mach number A^^ < 1 
and potential <I>a , and a downstream region characterized 
by Vb,A4b,^b- Here again, the wavelength of the pertur- 
bations is assumed to be longer than the size of the nozzle. 
The compact nozzle approximation enables us to compute 
the amplitude of the acoustic wave Sp~ /pa propagating 
upstream against the current, triggered by the refraction 
of an incoming acoustic wave 6p'^ /pa or by the advection 
of an entropy wave SSa- 

The following estimates of |7?,i | , | Qi | are deduced from Ap- 
pendix B: 
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(10) 



(11) 



1 -Ml Ci C^^n +VaCa' 

Following the argument of Appendix B, the most efficient 
entropic-acoustic coupling is expected at frequencies close 
to the refraction cut-off uj^-ut- The maximum value of \Qi\ 
is then of the order of 



|Q,|^0.3^^ ifc. 

Cl 



> Cl. 



(12) 



According to the Bernoulli Eq. (|^) with a gravitational 
potential: 

7- 1 /GAI 



^ = 2- 

son 111 

7+1 V ^s. 
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(13) 



The coupling coefficient \Qi\ deduced from Eqs. ( p^ ) and 
(O) can be very large if the sonic radius rson is close to the 
accretor, as in the radial accretion of a gas with an adia- 
batic index 7 close to 5/3 (Bondi 1952, Foglizzo & Ruffert 
1997). We wish to emphasize the fact that Eq. ( |l2[ ) comes 
out of a rather idealized toy model where both the refrac- 
tion of acoustic waves and their coupling to entropy per- 
turbations occur in a localized "compact" region. Never- 
theless, we assume it is realistic enough to expect | Qi | ^ 1 

if Cson > Cl. 

3. Entropy perturbation produced by a sound 
wave reaching a shock 

The acoustic energy propagating against the flow, released 
by the processes studied above, would escape far from the 
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accretor if no particular boundary is met. This energy 
however may feed a cycle if some physical process enables 
the conversion of acoustic energy into entropy perturba- 
tions. A similar cycle occurs in jet nozzles, where sound 
waves reaching the combustion region perturb the com- 
bustion and produce entropy perturbations. Here we focus 
on the case of a sationary shock satisfying the Rankine- 
Hugoniot conditions. For the sake of simplicity, we restrict 
ourselves to the case of a shock in the plane y, z with in- 
cident Mach number M.q in the x direction. The mach 
number M.s\y_ < \ immediately after the shock is related 
to A^o > 1 through the classical relation: 



Mi^ 



2 _2+(7-l)A^2 ^ ^_i 



"liMl -7 + 1 



> 



27 



(14) 



Let a sound wave Jp^J^ propagate against the stream along 
the X-axis in the subsonic region where the sound speed 
is Csh- The more general case of any angle of incidence 
is treated in Foglizzo (2000). The incident sound wave 
perturbs the shock surface and generates both a reflected 
sound wave (5p^ and an entropy perturbation 8S. We de- 
fine the reflexion coefficient 7?,sh and the efficiency of en- 
tropy production Qgh as follows: 



7^sh = 



Qsh 






(15) 
(16) 



A simple calculation in Appendix C, similar to Landau 
& Lifshitz (1987, chap. 90), shows that 7?.sh and Qsh are 
real functions of A^sh and 7. For a strong shock (A^g ^ 
2/(7 — 1)), the coefficients 7?.sh, Qsh depend only on the 
adiabatic index 7: 



7^sh 
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4. Global instability 

At a given frequency u, let us define t+, te as respectively 
the time of acoustic propagation and the advection time 
from the shock to the refraction point, and r_ the time 
of acoustic propagation from the refraction poin to the 
shock, taa = T~ 
cycle, while tea 
acoustic cycle. 



-r+ is the duration of the purely acoustic 
E t^^te is the duration of the entropic- 



J^.l. Instability of the entropic- acoustic cycle 

A first estimate of the entropic-acoustic instability is ob- 
tained by neglecting the effect of the purely acoustic cy- 
cle. The adiabatic advection of an entropy perturbation 





Fig. 2. Schematic trajectories of acoustic (dotted fines) 
and entropic (full lines) perturbations forming the 
entropic-acoustic and the purely acoustic parts of the cy- 
cle. Acoutic waves are reflected against the shock and 
refracted by the flow gradients with the respective effi- 
ciencies TZsh and TZa- Ta > rson is the refraction radius. 
Entropy perturbations are created by the interaction of 
acoustic waves with the shock with the efficiency Qsh. 
Their advection triggers outgoing acoustic waves with the 
efficiency Qa- Acoustic propagation times t-,t^ and en- 
tropic advection time te are also illustrated. 



of amplitude SS into regions of higher enthalpy produces 
acoustic waves Sp~^/psh estimated in Sect. g. When these 
acoustic waves reach the shock, new entropy perturbations 
are created and were estimated in Sect. g. This leads us 
to define the global efficiency Q of the entropic-acoustic 
cycle as follows: 



Q = QiQsh, 



(19) 



As for classical problems of modes in e.g. acoustic cavi- 
ties, the integral phase condition for the existence of an 
entropic-acoustic mode constrains the real part ujr,n of its 
frequency: 



riTT 
Tea ' 



(20) 



where the integer n must be even if Q > (e.g. in an 
accretion flow), and odd otherwise (e.g. in a jet nozzle). 
Note that the duration tea of the entropic-acoustic cy- 
cle also depends on the frequency through the position 
of the refraction point. The condition of instability of the 
entropic-acoustic cycle is |Q| > 1. This leads us to define 
the global growth (or damping) rate uii^n associated to the 
mode n: 



1 

Tea 



log|Q|. 



(21) 



4-2. Which accretion flows are unstable ? 



On the one hand, we know from Sect. 2.3 that the co- 
efficient jQij » 1 if Cson ^ Csh- On the other hand, we 
show in Appendix D that the coefficient Qsh rnay damp 
the instability if the shock is too weak. Thus the global 
efficiency |Q| exceeds unity for fiows such that Cson ^ Csh 
with a strong shock, as in the Bondi-Hoyle accretion with 
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Fig. 3. Growth rate as a function of the global efficiency 
Q, for taa/tea = 2/3, \R\ = 0.5. The dotted Hne corre- 
sponds to the entropic-acoustic cycle alone 7^ = 0. The full 
curves give the upper and lower bounds u^ of the growth 
rate when the contribution of the purely acoustic cycle is 
taken into account, w^ = for |Q| = 1 =F |7?.|. Eigenmodes 
are damped for uji < 0, and unstable for LOi > 0. 



The phase relation between the purely acoustic and 
the entropic-acoustic cycles is rather complicated (see 
Eq. E.5). The calculation of the eigenfrequencies requires 
to know the exact dependence of Q,TZ,tea^taa on the 
frequency. We show in Appendix E that the growth rate 
LOi is bounded to the interval [w" , uf] satisfying the equa- 
tion: 



\Q\e-'^' ^^■^±|7^|e- 



= 1. 



(26) 



The solutions uif of this equation are plotted as functions 
of Q in Fig. |. All the modes are unstable for I Q| > l-l-|7?.|, 
while all of them are stable for \Q\ < 1 — |7?.|. In the range 
1 — |7?.|<|Q|<1 + |7?.|, both stable and unstable modes 
may be found, depending on the phase of the acoustic 
cycle. A first order calculation of the contribution of the 
acoustic cycle when \TZ\/\Q\''^^^'^''^ < 1 leads to: 
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where n is even if Q > 0, and odd otherwise. The acous- 
tic cycle can be stabilizing or destabilizing depending on 
the phase miTAA/TEA- This phase dependence results in a 
wiggling in the spectrum of the eigenmodes, bounded by 



7 close to 5/3 at high Mach number. Following the com- 
pact approximation, the frequency of the most unstable 
mode is close to the refraction cut-off Wcut, which corre- 
sponds to a very high order rir, 

^cntTEA 



>1. 



(22) 



Eq. ( pl| ) indicates that the growth rate depends logarith- 
mically on I Q|, so that the maximum growth rate is always 
of the order of the inverse of the duration tea of a whole 
entropic-acoustic cycle if |Q| » 1. Using Eq. (p^), 

Wi.max ^ log -22!i, (23) 

Tea Csh 

4-3. Contribution of the purely acoustic cycle 

A refined estimate of the growth rate and frequency of 
the instability should take into account the sound waves 
^vth produced at the shock with the efficiency 7?.shj and 
their refraction near the sonic radius with the efficiency 
TZi. The global efficiency TZ of the purely acoustic cycle 
is: 



n = TliTZsh. 



(24) 



The eigenfrequencies must satisfy the following equation, 
deduced from Eqs. ( B.35 ) and (15-16|): 

Qg»u;rBA_^^git^r^,l ^ -^^ (25) 



4.4- Comparison with the particular case of jet nozzles 

The velocity and enthalpy profiles are quite different in a 
jet nozzle and in an accretion flow. The gas is accelerated 
in a nozzle at the expense of the enthalpy, whereas gas 
flowing into a gravitational well <1> = —GM/r is both ac- 
celerated and heated. The limited variation of enthalpy in 
a nozzle, deduced from the Bernoulli equation with <f> = 0, 



,2 
'son 
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2 + {j-l)Ml 
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<1, 



(29) 



results in a moderate efficiency of the extraction of acous- 
tic energy from the advection of entropy perturbations 
(|Q| < 1). Nevertheless, the combustion zone plays a 
much more active role than the shock described above, 
since the energy input due to the chemistry of the com- 
bustion enables the actual efficiency |Q| to exceed unity. 
In this respect, the "rumble" instability of ramjets (Abou- 
seif, Keklak & Toong 1984) is different from the entropic- 
acoustic instability described here. Both are based, how- 
ever, on the cycle of entropic and acoustic perturbations in 
the subsonic region of the converging flow. In an accretion 
flow, the efflciency |Q| increases with the frequency up 
to the refraction cut-off, thus favouring a high frequency 



instability 



:$> 1). By contrast, the most unstable 



modes in Abouseif, Keklak & Toong (1984) are the low 
frequency modes. 



T. Foglizzo & M. Tagger: Entropic-Acoustic instability 



5. Discussion 

We have described how the interplay of entropic and 
acoustic waves in a subsonic flow between a shock and 
a sonic point might lead to an entropic-acoustic linear in- 
stability. The ingredients for this type of instability are 

(1) an accelerated subsonic flow with a significant in- 
crease of sound speed (i.e. of temperature), 

(2) an outer boundary, such as a shock, capable of con- 
verting acoustic waves into entropy perturbations. 

The stability criterion is directly related to the calcula- 
tion of a global parameter Q describing the coupling of 
entropy and acoustic waves in the flow. The effect of the 
purely acoustic cycle is negligible if |Q| S> 1, but becomes 
important near marginal stability (Eq. ^). Although a 
detailed stability calculation is required for any particular 
flow topology, general arguments based on the compact 
nozzle approximation suggest that the entropic-acoustic 
cycle is potentially unstable if the sound speed at the sonic 
point significantly exceeds the sound speed immediately 
after the shock (Eq. g3|). The two timescales involved in 
the entropic-acoustic instability are: 

- the period of the most unstable pertubations, of the 
order of the refraction cut-off, 

- the growth time, comparable to the time needed to 
advect an entropy perturbation from the shock to the 
sonic radius. 



The entropic-acoustic cycle could then occur between this 
shock and the inner sonic surface. On a larger scale, the 
gas falling from the Lagrangian point Li onto the external 
radius of an accretion disc might also produce a shock ca- 
pable of closing the entropic-acoustic cycle. These possible 
applications deserve a careful analysis which is beyond the 
scope of the present paper. 

5.2. Towards a 3-D description of the instability 

Our one-dimensional description assumed that all the out- 
going acoustic waves reach the shock, and that all the 
induced perturbations of entropy reach the sonic surface. 
Depending on the topology of the flow, possible losses from 
the entropic-acoustic cycle must be taken into account. In 
the case of Bondi-Hoyle accretion, outgoing acoustic waves 
reaching the shock out of the accretion cylinder produce 
entropy perturbations which will not be accreted, thus di- 
minishing the effective efficiency Q of the cycle. A 3-D 
description of the entropic-acoustic cycle should also take 
into account the effect of vorticity perturbations, which 
are naturally created when an acoustic wave perturbs a 
shock with a non vanishing angle of incidence. Like the 
entropy perturbations, the advection of a vorticity per- 
turbation in a converging flow excites acoustic waves (e.g. 
Howe 1975). Their destabilizing effect will be studied in 
the particular case of radial accretion (Foglizzo 2000). 



5.1. Astrophysical applications 

The entropic-acoustic instability might be at work in sev- 
eral astrophysical situations: as already mentionned in the 
introduction, the supersonic motion of a compact star in 
a gas leads to the formation of a bow shock ahead of the 
accretor. The development of the entropic-acoustic insta- 
bility in the subsonic region between the shock and the 
sonic surface is very likely, since Cson/csh ^ 1 in these 
flows (e.g. Ruffert 1996 and references therein). A detailed 
analysis will be presented in a forthcoming paper. 
The acoustic-entropic instability might also be present in 
some accretion discs. Most stability calculations concern- 
ing ADAFs focussed on the local thermal stability of the 
disc (Narayan & Yi 1995, Abramowicz et al. 1995, Kato et 
al. 1996, 1997, Wu & Li 1996). By contrast, the entropic- 
acoustic instability is a global instability which depends 
on the outer boundary condition of the flow. The excita- 
tion of acoustic waves from an advected entropy perturba- 
tion in an ADAF has already been observed in numerical 
simulations by Manmoto et al. (1996), which include the 
additional effects of radiative cooling (bremmstrahlung) 
and viscous heating. Their work was concerned with the 
time delay problem of Cygnus X-1 and did not involve an 
outer boundary condition capable of closing the entropic- 
acoustic cycle. Some global disc solutions match a thin disc 
to an inner ADAF through a shock (Lu, Gu & Yuan 1999). 



5.3. Observable consequences 

Understanding the fundamental mechanism of linear in- 
stability is only the very first step towards a possible iden- 
tification of this mechanism in the astrophysical observa- 
tions. According to our linear description of the entropic- 
acoustic instability, we expect to observe significant vari- 
ations of the mass accretion rate and displacements of 
the shock surface as the acoustic energy increases between 
the shock and the sonic surface. The determination of the 
spectrum of eigenfrequencies should give valuable infor- 
mations about the timescale of the most unstable modes. 
Numerical simulations should be the most direct way to 
describe the ultimate non-linear development of this insta- 
bility. We saw that a correct description of the instability 
requires an accurate propagation of acoustic and entropy 
waves between the shock and the sonic radius. This could 
be a serious numerical difficulty for flows where Cson ^ Cgh, 
since the frequency of the most unstable perturbations 
could exceed the growth rate by several orders of magni- 
tude. The numerical simulations of Bondi-Hoyle- Lyttleton 
accretion (e.g. Ruffert 1995) suggest that this instability 
might lead to some kind of turbulence. This could also be 
particularly promising if the same mechanism also applies 
to accretion discs. 
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Appendix A: Energy release of a localized 

perturbation of entropy through a 
nozzle 

We consider in the region I an entropy perturbation SS, cor- 
responding to a perturbation of mass 5m initially in pressure 
equilibrium with the surrounding gas {5pi — 0, 5vi — 0). The 
time scale r associated to this perturbation of entropy is re- 
lated to its length I: 



I 

1 

V 

Sp 



5m 



Mot, 



where Mo is the unperturbed mass flux. The entropy pertur- 
bation is conserved when advected: 



5S : 



■ log 1 



5p 
P 



5p 
P 



(A.3) 



The advected perturbat ions of density and sound speed {(? = 
"ip/p) follow from Eqs. (A.3): 



5pi 
5cf 



= exp 



= exp 



7-1 
7 
7-1 



5S 



7 



5S - 1. 



(A.4) 
(A.5) 



The energy flux Fq + 5F-i in the perturbation of entropy is 
defined according to La ndau & Li fshit z (1987, Chap. 6), and 
deduced from Eqs. (§), (^) and (^): 



5Fi = {Mo + 5M){Bo + 5B)-MoBo, 



= Mo I Y + *i 



exp 



7-1. 

7 



5S -1 



(A.6) 
(A.7) 



Replacing the index 1 by the index 2 gives the energy flux Fo + 
5F2 in the entropy perturbation once it has reached pressure 
equilibrium with the region II. The total energy 5Ei2 released 
during the advection through the nozzle is therefore: 



5Ei2 = {5Fi-5F2)t, 
= Mor{h2 - hi) 
= (h,2 — hi)5m. 



exp 



7-1. 



5S -1 



(A.8) 

(A.9) 

(A.IO) 



Appendix B: The compact nozzle in a potential 

*(r) 

B.l. Conservation equations 

The variations 5M and 5B of the mass flux M and Bernoulli 
parameter B within a pressure perturbation 5p^ (±1 depend- 
ing on the direction of propagation of the perturbation) is to 
first order: 



5M^ 
M 



5v^ M±l 
7P M 



5B^ = c\l±M)^-^ 
IP 



(B.l) 
(B.2) 



In an entropy perturbation 5S in pressure equilibrium with its 
surrounding, the corresponding variations Shf^ and SB"^ are as 
follows: 



5M^ 
M 

5B'' 



7-1 



5S, 



-5S. 



7 



(B.3) 



(B.4) 



Writing the conservation of the perturbed entropy, mass fiux 
and Bernoulli parameter across a subsonic nozzle leads to the 

(B.5) 
(B.6) 
(B.7) 

The case of a supercritical nozzle requires one more equation 
describing the effect of perturbations on the sonic surface. The 
Bernoulli equation can be written at the sonic point rson in a 
fiow tube of entropy S with mass flux M and local cross section 
As-. 



(A.l) 


equations for the subcritical nozzle: 


(A.2) 


5Sa = 5Sb, 

5M+ + 5M- + 5M: = 5M+ + 5M- + 5 Ml, 


srtur- 


5Bt + SB- + 5BI = 5B+ + 55" + 5Bl 



B - ^(r,, 



77-1-1 
27-1 



Me^ 



2(t-1) 



72A 



(B.i 



Differentiating the Bernoulli equation and the mass conserva- 
tion at the sonic point leads to 

5A 



■son -7- = 5$(rs, 



(B.9) 



Combining this equation with the differential of Eq. (B.8) pro- 
vides us with the additional equation required to solve the case 
of a supercritical nozzle: 



5M ,^ 5B „ 

M Clon 



(B.IO) 



B.2. Incident acoustic wave propagating with the stream 
Subcritical nozzle: 

5pZ _ 1 + Ma VbCb - VgCg 5p+ 
7Pa 1 - Ma VbCb + VaCa ^Pa ' 

Spt l + Ma 2Mbcl 5pt 



'yPb 1 + Mb VbCb + VaCa ^Pa ' 

Supercritical nozzle: 

Sp~ l+Ma cLn - VaCa 5p+ 



(B.ll) 
(B.12) 



IPa 


1 - Ma clan + VaCa -JPa ' 


Spt 


1 -1- Ma cl Cgon -1- VbCb 5p+ 


IPb 


1 -1- Mb cl don + VaCa ^Pa 


5Pb 


1 -1- Ma cl Cson " VbCb 5pt 



'yPb I - Mb cl Cion + VaCa ^Pa 



(B.13) 
(B.14) 
(B.15) 



B.3. Incident entropy perturbation advected with the 
.stream 



Subcritical nozzle: 



Ma 



5S 



5Pa ^ 

JPa 1 - Ma VbCb + VaCa 7 

Spt ^ Mb cl - d 5S 

JPb 1 + Mb VbCb + VaCa 7 



(B.16) 
(B.17) 
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Supercritical nozzle: 
Ma cl 



IPa 

&pb 

ipb 



5S 



7 



(B.18) 



■^\Cson H~ VaC-a) 



2(cion + VaCa) 






< 1+7W„ c^ 



5S_ 

1 

7 
7 ' 



(B.19) 



(B.20) 



B.^. Acoustic flux F^ propagating upstream 

Since our calculation is based on a linear approximation, we 
may define efficiencies 7?.a , Qa of the compact nozzle as follows 



SPa Tp Spt 

Pa Pa 



(B.21) 



where TZ a, Qa are ea sily d educed from Eqs. (B.ll) and (B.16) 
or Eqs. (B.13) and (B.16), depending on whether the nozzle 
is subcritical (Mb < 1), or supercritical {Mb > !)• The pres- 
sure perturbation Spa g enerated u pstream in the case of a 
supercrit ical nozzle (Eqs. BT3-B.18) corresponds to the limit 
of Eqs. (B.11-B.16) when cj, = Csan, the sound speed at the 
sonic point. This is consistent with the fact that no information 
about the fluid properties beyond the sonic p oint m ay reach 
the subsonic region upstream. Note from Eq. (B.18) that the 



sign of Qa depends on whether the sound speed of the acceler- 
ated gas increases or decreases. The Bernoulli equation implies 
that Qa < in a classical nozzle ($ = 0), while Qa > in 
an accretion flow. The flux F^T of acoustic energy propagating 
against the stream when an acoustic wave carrying the acous- 
tic flux Fa (with 5S = 0) prop agate s into a supercritial nozzle 
is deduced from Eqs. (|6|) and ( B.13 ): 



ni 



Fa 

F^ 



- VaCa 



Cson V VaCa 



(B.22) 



where the refraction efficiency TZa is defined i n ter ms of the 
acoustic energy, while TZa was defined in Eq. (B.21) in terms 



of the amplitude of pressure perturbations. It is remarkable 
that the refraction Fa I Fa is independent of the frequency of 
the incoming acoustic wave. This is true only at low frequency, 
as long as the compact approximation is satisfied. Refraction 
is expected to become negligible (Fa I Fa ^ 1) at high fre- 
quency, when the wavelength of an incoming acoustic wave is 
much shorter than the lengthscale A of the flow inhomogeneities 
(i.e. the nozzle size). The refraction cut-off ijjcut is deflned as 
the maximum frequency leading to a significant refraction (e.g. 
FalFi ^\I2). 

The flux Fa of acoustic energy propagating against the stream 
when an entropy perturbation &S (with &p^ = Oj is ad yecte d 
into a supercritical nozzle is deduced from Eqs. dq) and (B.lq): 



Fa = Ml^VaCa 



+ VaCa 



(B.23) 



Here again, the acoustic flux Fa is independent of the fre- 
quency of the incoming entropy wave, as long as the compact 
approximation is satisfied. Numerical calculations by Marble & 



Candle (1977) show that the efficiency of sound emission from 
an entropy perturbation also decrea ses to zero if its wavelength 
is much shorter than A. Eq. (B.23) is rewritten using the re- 
fraction coefficient TZa'- 

2 / „„\ 2 



Fa 



— MoCsonTla 



l~TZa 



The ratio (c,, 



i ~T l\'a / \ Cson CaVa 

- Ca)/(cson - CaVa) is bouuded by 

1 for Ca < Cson, 

^^'sonCson 



— •(B-24) 



24, 



for Ca ~ Cs, 



(B.25) 
(B.26) 



where A^son, Cson are the derivatives of the mach number and 
sound speed at the sonic point. Schematically, TZa is constant 
for cij <^ tJcut and decreases to zero for ijj ^ oicut. The maximal 
value of 7^a(l-7^a)/(l-l-■R.a) ~ 0.09 is reached for 7?.a ~ 0.6, i.e. 
naCa '~ Cgon/4. We couclude that the acoustic flux propagating 



against the stream when F2 

2 



is bounded by: 



Fa < 0.09 Moci 



(B.27) 



B.5. Combination of a slowly accelerated flow and a 
compact nozzle 

Let us build a toy model of a flow Mi < M < 1, especially de- 
signed to allow a simple calculation of the interplay of acoustic 
and entropic perturbations of frequency to. This flow is com- 
posed of: 

- a region of propagation Mi < M < Ma where both 
the entropic-acoustic coupling and the acoustic refraction at 
the frequency ui are negligible, and where the flow is slowly 
accelerated, 

- a compact region Ma < M < 1 of coupling, approxi- 
mated as a supercritical compact nozzle. 

The index a notes the entrance of the compact nozzle. These 
two approximations are compatible if the lengthscale of the 
flow inhomogeneities is small compared to the wavelength of 
the perturbations in the compact region, and large compared 
to this wavelength in the propagation region. The real coeffi- 
cients Pi , T'e describe the amplitude variation of the pressure 
and entropy perturbations during their propagation between 
Ml and Ma- 

SpI ,^r^ 

pi 

(5n+ ■ 
e 

Pa 

VESSie"^^". (B.30) 

The values of T-'± are directly deduced from the conservation 
of the acoustic flux from Ma to Mi (F^ ~ Fa), and T'e is 
deduced from the conservation of entropy: 



5pi_ 

Pa 
SPi 
Pi 
SSa 



V+ 



V- 



(B.28) 



(B.29) 



r+ = 



P- 



1+Mi 


MlCi 


1+Ma 


MjCa 


1-Ma 


MjCa 



1-Mi 



Mici 



Ve 



(B.31) 

(B.32) 
(B.33) 
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Combining Eqs. (B.21) and (B.31-B.33), we obtain the ampli- 
tude Sp^ of the acoustic waves propagating against the stream 
at the position A^i, as a function of the incident perturbations 
5pf,5Si: 



7^l = V-V+TZae'^^^^ 



(B.34) 
(B.35) 




a re de duced from Eqs. (B.l^), (BJJ) and 



i.31| ) to (B.33). From the discussion of Sect, 
strongest acoustic flux F{' is obtained for TZa ~ 0.6 



B.4 



the 



Appendix C: Reflexion of a sound wave at a 
shock 

Consider an adiabatic shock with incident Mach number A^o 
in the x direction. Let a sound wave propagate against the 
stream in the subsonic region where the sound speed is Csh , with 
wavevector k- in the frame moving with the fluid. Its frequency 
oj in the frame stationary with respect to the shock obeys the 
following dispersion relation (e.g. Landau & Lifschitz, 1986, 
chap. 68): 

{uj - kvshf = cl^k^ . (C.l) 

The two roots fcip of this equation correspond to the incident 
and reflected acoustic waves. The incident acoustic wave per- 
turbs the shock surface and generates both a reflected acoustic 
wave and an entropy perturbation. The pressure perturbation 
5ps\i and the entropy perturbation 5S axe obtained by writing 
the jump conditions for an adiabatic shock perturbed with the 
velocity A« in the x direction: 

47-Mo 



&PBh 

Psh 

5S = - 



Av 

2'yMl -7-1-1^' 
2[l+ Ml) Av 



(7- 



l)Ml vo ' 

A-fiMl - 



Av 



{2 + {'y-l)Ml){2jM'o-l + l) vo 



(C.2) 
(C.3) 

(C.4) 



The perturbation is decomposed onto the ingoing and outgoing 
acoustic waves, and the advected entropy perturbation 



SpBh 
Svsh 



SPsh + Spth^ 
Sv' + SvX + Svti 



(C.5) 
(C.6) 



In addition to the dispersion equation (C.l) for acoustic waves. 



the acoustic and entropy perturbations are described by: 



(a.- 



k±)Svi 


7 Psh 


LU 


= ke Vs\\ 5 


kjvth 


= 0. 



After some algebra with Eqs. (0.1) to (C.9), we obtain: 
2 3 - 7 - 2(7 - 1)Ma, 



'T^sh — 



Qsh 



Mb 



(1 



1 + A4sh 

- TWeh)" 



7 + 2(7- 
4 



l)Xsh 

>0 



<0, 



7^.l 



X,h 3-7-^2(7- 1)X 



l-.Msh 1 3 - 7 - 2(7 - 1)Mby. 
l + X,hJ 3 - 7 + 2(7 - l)7Meh 



(C.7) 

(C.8) 
(C.9) 

(CIO) 
(C.ll) 

(C.12) 



where V^^^ is the reflection coefficient of acoustic waves defined 
in terms of energy fiux rather than amplitudes. For a strong 
shock, iZsYi is an increasing function of 7 {iZs\i ~ —0.14 for 
7 = 5/3): 



22 --y2(7-l)! 

22 — ^2{j — 1)1 



(C.13) 



Appendix D: Estimates of the global efficiencies 
Q and Ti. based on the compact 
nozzle approximation 



7^ = V^V+TZaTlBh, 

— TZalZsh- 



(D.l) 
(D.2) 



TZ depends on the frequency through TZa'- TZ ~ 7?.sh for cu ^ 
tJcut , and 7?. ~ for OJ » Ucut . 



Q = P-pEQaQsh, 



/ VgCg \ 

VWshCsh/ 



+ VaCa 



a f 
/ah, 



.fsh = 



-^sl 



Xsh 



3 - 7 + 2(7 - l)>1,h 



(D.3) 
(D.4) 



(D.5) 



The efficiency Q decreases to zero as (1 — A^sh) if the shock 
is weak. For a strong shock, /sh is a decreasing function of 7 

(/,h ~ 1.7 for 7 = 5/3): 



/sh 



27 



22 



7-ly 2t-H7t(7-l)t 
Q is rewritten using the refraction coefficient TZa'- 



Csh 



1-TZa 
1+TZa 



■ CaVa 



/si 



(D.6) 



(D.7) 



Since Qa depends strongly on the frequency, the efficiency Q 
is a function of uj. Schematically, iZa ~ 1 for uo <^ cjcut, and 
TZa ~ for u) ^ Wcut, so that Q reaches its maximum value 
for a frequency close to the cut-off frequency. In particular 

Qmax > 1 if Cson > Csh: 



son "--son 



/sh > ymax > U.J/sh t H TTT 

Csh \ 2c, 



^:^. (D.l 

Csh 



Appendix E: Contribution of the acoustic cycle 

The complex frequency to is decomposed into its real part Ur 
and imaginary part uui. Let us introduce the parameter a; > 
and transform the eigenmodes equation as follows: 



Q_ 

\Q\ 



1 

Tea 
1-1 



+ X 



loglQil, 



\Q\ 



^EA 



(E.l) 
(E.2) 



The solutions of Eq. (E.2) are bounded by the interval X- < 
X < x+, where X- < 1 < x+ axe solutions of the equation: 



1-^ |7^l 



l±x^ "-^-^ 



\Q\ 



^AA 
''EA 



= X±, 



(E.3) 
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from which we obtain Eq. (Eq) satisfied by t he cor resp onding 
growth rates ujr^. We deduce from Eqs. (E.l) and (|E.3|) that 



^nlQ| = lT|7^|)=0. (E.4) 

The real part u;,. must satisfy the following phase equation, 



deduced from Eq. (E.2) 



SmUJr{TAA —TEA 



Sm UJrTAA 



-- ^^SmUr{TEA~rAA)^ (E.5) 

SaiLUrTEA 

^^QsinU^riTAA-TEA)^ (E.6) 

\Q\ SmLJrTAA 

If the parameter \Ti\/\Q\'^^^''^^^ <^ 1, we can estimate the first 
order contribution of the purely acoustic cycle to the entropic- 
acoustic instability: 
7^ 






cos LUrTAA- (E-7) 



Together with Eq. pHj - lE^ , we obtain Eqs. (0-||). 
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